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Abstract 



The UHF wave function may be written as a 
spin-contaminated pair wave function of the 
APSG form, and the overlap of the alpha and 
beta corresponding orbitals of the UHF solu- 
tion can be taken as a proxy for the strength 
of the correlation captured by breaking sym- 
metry. We demonstrate this with calculations 
on one- and two-dimensional hydrogen clusters 
and make contact with the well studied Hub- 
bard model. The UHF corresponding orbitals 
pair in a manner that allows a smooth evolution 
from doubly occupied orbitals at small distance 
to one in which wave function breaks symme- 
try, segregating the a and /5 electrons onto dis- 
tinct sublattices at large distances. By perform- 
ing spin projection on these UHF solutions, we 
address strong correlations that are difficult to 
capture at intermediate distances using a sin- 
gle determinant. Approved for public release: 
LA-UR-13-22691. 

1 Introduction 

The proper treatment of what has come to be 
known as the "strong correlation" problem re- 
mains a fundamental challenge to our under- 
standing of unconventional superconductivity, 
frustrated spin lattices, the heavy fermion prob- 
lem, and the class of materials known as Mott 
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insulators. Materials are typically considered 
strongly correlated when conventional density 
functional approximations fail to qualitatively 
describe their properties. For example, the lo- 
cal density and generalized gradient approxi- 
mations (the LDA and GGA) of density func- 
tional theory (DFT) predicted NiO to be a fer- 
romagnetic metal. This catastrophically failed 
to describe the experimental properties: NiO 
behaves as an antiferromagnetic insulator with 
an optical gap of the order of 4 eV. ^ 

This led to statements by some that "DFT 
cannot describe Mott insulators." Recent re- 
search has shown that this failure lies in the lo- 
cal or semi-local exchange-correlation approx- 
imations typically employed in the field, and 
not with DFT itself. Including a fraction of 
the fully non-local HF exchange term in the 
functional remedies many of these problems, 
and these hybrid functional yield significantly 
improved band gaps, lattice constants, and 
magnetic properties for Mott insulators.'^ The 
strong correlations responsible for opening a 
gap in NiO do not vanish if it is made metallic 
by application of pressure, and so conventional 
local and semi-local approximations to DFT 
have problems with strongly correlated metals 
as well. Mounting evidence suggests that the 
hybrid functionals do not perform well in this 
regime either .I^EEI 

Theorists have struggled with this problem for 
decades, and many approaches that build upon 
traditional DFT have attempted to address 
it. These include the self-interaction-correction 
(SIC) to DFT,^ DFT+U approaches ,133 and 



many-body approaches such as the GW ap- 
proximation,'^*^ and dynamic mean-field the- 
ory (DMFT).^ The SIC offered a signifi- 
cant improvement over conventional DFT in 
the Mott insulator regime but vanishes for a 
metal, necessitating additional approximations 
in that realm. ^^ The GW approach, a many- 
body Green's function method based on the 
random phase approximation, successfully pre- 
dicted band gaps in many materials and some 
Mott insulators. However, the strength of this 
method, particularly its focus on the quasipar- 
ticle excitation spectrum, also made it difficult 
to obtain a total energy even though this energy 
is connected to the RPA ground state correla- 
tion.!^ DMFT, while formally a non-empirical 
embedding method, relies in practice on mate- 
rial dependent parameters for its utilization, as 
do the DFT+U approaches .l^l A truly predictive, 
parameter free, ab initio description of the elec- 
tronic structure of strongly correlated materials 
remains a significant challenge for theory. 

The Mott transition bears striking similar- 
ities to one of the oldest problems in quan- 
tum chemistry: making and breaking chemical 
bonds. Specifically, consider breaking the sim- 
plest bond, that of H2 (see Figure [l]). Near 
equilibrium, a symmetry restricted Hartree- 
Fock (RHF), doubly occupied ag molecular or- 
bital (MO), dominates the wave function. This 
description fails abysmally in the dissociation 
limit; the double occupancy implicit in the 
a'^ configuration contributes energetically un- 
favorable ionic components to the wave func- 
tion at large bond lengths. In this limit, the 
bonding and anti-bonding molecular orbitals 
become nearly degenerate; when expressed in 
this symmetry restricted MO basis, a proper 
wave function requires two configurations, a^ 
and al with approximately equal contributions 
(the full configuration interaction, or FCI wave 
function in a minimum basis). A superposi- 
tion of these two configurations with variation- 
ally determined weights evolves smoothly from 
the delocalized MO limit to a localized valence 
bond (VB) description. Similarly, the tran- 
sition from a simple metal to a typical Mott 
insulator requires a wave function which con- 
tinuously evolves from the independent elec- 
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Figure 1: Dissociation curve for H2 computed 
with the ST0-3G basis set for illustration pur- 
poses. The RHF single determinant qualita- 
tively describes the molecule near equilibrium 
but does not dissociate correctly (see text). By 
artificially breaking spatial and spin symmetry, 
the single determinant UHF solution yields a 
qualitatively correct curve in both limits. 

tron, delocalized Bloch orbital description into 
a strongly correlated, valence bond-like wave 
function as the localization increases. 

The success of hybrid functionals in describ- 
ing Mott insulators relies, in part, on relaxing 
symmetry constraints imposed on the molecu- 
lar orbitals. In the unrestricted Hartree-Fock 
(UHF) or Kohn-Sham (UKS) approximation 
each electron occupies its own spatial orbital.!^ 
This allows the wave function to localize the 
electrons when appropriate, one to each atom 
in H2 at large distance or one per site in an 
anti-ferromagnet, thereby avoiding the unfavor- 
able ionic component implicit in the restricted 
approximation. Relaxing symmetry constraints 
sacrifices the quantum numbers associated with 
both spatial and spin symmetries. Neverthe- 
less, the unrestricted (UHF) wave function al- 
lows a qualitatively (if not quantitatively) ac- 
ceptable description of both the Mott insula- 
tor and simple metal limits. We should not 
be surprised that hybrid functionals have diffi- 
culty describing the correlated metal regime,'^'^^ 
corresponding to the intermediate bond lengths 
of Figure [1} In this region, strong correlations 
expose the problems inherent in single deter- 
minant methods, and a successful theory must 
delicately balance the competition between lo- 
calization and delocalization. 



In this work, we focus on understanding this 
intermediate "recouphng" region, in which nei- 
ther the MO nor VB picture, neither the re- 
stricted or unrestricted approximations are ap- 
propriate. We begin by reviewing some con- 
cepts regarding broken symmetry and pair wave 
functions, then briefly discuss wave functions 
which cannot be expressed in terms of pair func- 
tions. Next we describe projection methods 
which restore the symmetries, iUust rated with 
calculations on flnite size hydrogen networks. 
Finally, we conclude by discussing the implica- 
tions for truly periodic networks. 

2 Pair Wave Functions 

We begin by reviewing the wave function of 
diatomic hydrogen as a function of bond dis- 
tance. Our discussion utilizes a minimal ba- 
sis set, the simplest expansion which captures 
the near degeneracies between the molecular 
orbitals that occur as the bond is stretched. 
A minimum basis set also provides a concep- 
tual link with the Hilbert space utilized in the 
widely studied Hubbard model,'^ though we re- 
tain the exact Hamiltonian in the calculations 
to follow. These near degeneracies give rise to 
our primary concern, the so-called left-right, 
or static correlation. Static correlation influ- 
ences the form of the wave function, and can 
generally be accounted for in the Hilbert space 
of the minimum basis. In a condensed mat- 
ter context this might be referred to as antifer- 
romagnetic correlation. Dynamic correlation, 
or the instantaneous avoidance of one electron 
by another, constitutes a much weaker effect 
and, generally speaking, the local and semi- 
local correlation functionals of DFT describe 
it adequately. That is to say, for a molecule 
qualitatively well described by a single deter- 
minant of doubly occupied orbitals then many- 
body perturbation theory, coupled-cluster the- 
ory, etc., can describe the dynamic correlation, 
although it requires an extension of the basis 
set. 

To this end, consider a single Is orbital cen- 
tered on each of two hydrogen atoms, denot- 
ing these as 0/ and 0^, assigned to the left and 



right atomic centers, respectively. For simplic- 
ity, we assume orthogonality of these atomic 
wave functions, though this is not necessary 
in principle. Form molecular orbitals using the 
standard linear combinations of atomic orbitals, 
giving gerade (a^, bonding) and ungerade (a^^, 
anti-bonding) orbitals: 
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The RHF determinant with the bonding orbital 
doubly occupied approximates the exact wave 
function quite closely at the equilibrium bond 
length. This wave function can be written as: 

^[a,(l)«(l)a,(2)/3(2)]=^[a,a,], (2) 

where A is the anti-symmetrizer, a and /3 are 
the normal spin functions and we have used an 
overline to indicate beta spin-orbit als. One may 
then expand this wave function in the original 
atomic orbital basis: 

a,a, = |$J^)|0,0) + |$J)|0,0), (3) 

by defining the following spatial functions: 



\^l) = -7^ (Ml + (I^Av) 

and the singlet spin wave function: 
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These two spatial functions describe the two 
components of gerade symmetry that con- 
tribute to the ^Tig ground state of H2. The flrst 
term, |$^), describes the covalent portion of 
the wave function, in the Heitler-London (HL) 
sense. The ionic component, |$^), describes the 
symmetric combination (S^) of charge-transfer 
terms H^H~ and H~H^. For completeness, de- 
flne two spatial wave functions of ungerade sym- 
metry that play roles in the excited ^S^ and ^S^ 
states and appear later in the discussion of the 



UHF wave function: 
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At long range, the correct wave function for 
the ground state of H2 contains only the HL 
term, |$?). The RHF wave function dissociates 
improperly precisely because the ionic piece im- 
poses an energy cost equal to the ionization po- 
tential of one hydrogen atom minus the elec- 
tron affinity of the other. This quantity is re- 
ferred to as U in the Hubbard ubiquitious in 
the condensed matter literature and 7^^ in the 
Pariser-Parr-Pople model familiar to quantum 
chemists. It is a simple measure of the on-site 
electron-electron repulsion. 

Qualitatively, what does this mean? At small 
distances, delocalizing the electrons over both 
centers allows the electrons to gain sufficient 
kinetic energy to overcome the on-site repul- 
sion. At large distances, the kinetic energy ob- 
tained from delocalization becomes small, over- 
whelmed by the associated on-site electron re- 
pulsion, and the electrons localize. An accept- 
able approximation requires a wave function 
with can evolve smoothly between these two 
limits. 

2.1 The Generalized Valence 
Bond Wave Function 

The simplest wave function for H2 that 
smoothly evolves from a delocalized MO de- 
scription to a localized VB description is the 
generalized valence bond (GVB) wave func- 
tion.'^^'^For two electrons this may be thought 
of as a multi-determinant extension of RHF to 
an open-shell singlet pair (geminal) wave func- 
tion of the form: 



'^GVB = A[XlXr]\0,0), 



(7) 



where the non- orthogonal spatial functions Xi 
must be determined self-consistently at each 
distance. The non-orthogonality allows each or- 
bital to polarize continuously as a function of 
bond length, providing enough variational ffex- 



ibility to successfully break the bond. 

Expanding each of these non-orthogonal or- 
bit als in the original MO basis: 
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(8) 



introduces a mixing parameter A, a real number 
between -1 and 1. The spatial overlap of the 
non-orthogonal orbitals can be calculated: 



Sir = iXllXr) 



1-A 



1 + A' 
and after some algebra, Eq.Q becomes: 

^(1 + A2)' ^^' ' ^ 
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This form makes explicit the relationship be- 
tween the ionic and covalent portions of the 
wave function, while remaining equivalent to 
the multi-configuration representation more fa- 
miliar to quantum chemists: 
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[a^^-Aa^]|0,0). (11) 



For a single pair in a minimal basis set, this 
wave function is exact, with the mixing pa- 
rameter A and the orbitals determined self- 
consistently. 

We have belabored this simple two-electron 
example because we wish to emphasize the use- 
fulness of this construct for strongly correlated 
materials. Pushing the solid-state analogy fur- 
ther, this simple pair function continuously con- 
nects an ionic charge-density- wave (A = — 1), a 
simple metal (A = 0), and the antiferromag- 
netic singlet (A = 1) with a single mixing pa- 
rameter. 

To extend this pairing concept to systems 
with more electrons, one writes the wave func- 
tion for A^e electrons as an anti-symmetrized 
product of geminals (APG). For an extensive 
review of geminal wave functions see, for ex- 
ample, the work of Surjan'221 and the references 



therein. In the APG formahsm, electrons not 
exphcitly correlated are grouped in the core and 
frozen from excitation. Each correlated MO ipi 
pairs with an unoccupied partner ipif in the fol- 
lowing manner: 



N 

ip = A[core] 1 [ 

i=l 
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If one imposes the condition of strong orthog- 
onality between the pairs, the wave function 
is denoted as an anti-symmetrized product of 
strongly orthogonal geminals (APSG). In ad- 
dition to strong orthogonality, enforced by re- 
stricting each pair to distinct occupied and vir- 
tual orbital spaces, one may further restrict the 
pairs to pure spin states, a constraint denoted 
as "perfect pairing". ^^ Such a wave function ex- 
plicitly neglects inter-pair correlation effects, as 
well as same-spin correlation, that may be im- 
portant in many cases, particularly for high spin 
multiplets. 



2.2 Broken Symmetry 
Functions 



Wave 



The UHF approach eschews the multi- 
determinant nature of GVB in favor of a single 
determinant wave function. Allowing differ- 
ent spins to occupy different spatial orbitals 
provides enough flexibility for the electrons to 
balance the contributions from the covalent and 
ionic terms. At short distances this approach 
yields the RHF wave function, but as the bond 
stretches past some critical distance, known 
as the Coulson-Fischer (CF) point, the kinetic 
energy and nuclear attraction terms associated 
with the ionic portion no longer compensate 
for the on-site repulsion, and a distinct broken- 
symmetry UHF determinant emerges as the 
ground state. Artificially breaking spatial sym- 
metry by placing electrons of different spin in 
different orbitals finesses the electron correla- 
tions that dissociate the bond, projecting away 
the ionic component of the RHF wave function. 
Breaking symmetry in the presence of orbital 
near degeneracies should not lower the energy, 
but the fact that the UHF energy lies lower than 
the RHF energy lies at the heart of Lowdin's 



symmetry dilemma: one obtains a better ap- 
proximation to the total energy by artifactually 
sacrificing good quantum numbers. ^^ 

Returning to the simple two-electron case, the 
UHF wave function takes the form: 
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where xi ^^d Xr are the non-orthogonal orbitals 
defined by the mixing coefficient A in Equation 
([S]). In the dissociation limit, these orbitals re- 
duce to the atomic orbitals associated with each 
atom, and the wave function becomes either: 



or 



'^(Pr^ 



both of which break spatial symmetry. The spin 
products a/3 and I3a are not eigenfunctions of 
spin, but rather a linear combination of the sin- 
glet and triplet functions: 

^ |0,0) + |1,0) ^ ^^ , 

a/3 = ' ^ ^' ^ and (14) 
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where |1,0) is the m^ 
triplet: 

|i,o) 



(15) 



component of the 
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Rewriting the full UHF wave function by sep- 
arating the spatial parts into symmetric and 
anti-symmetric components, we identify the 
spatial functions associated with the singlet and 
triplet spin components: 
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K\hO). 



The singlet component (up to normalization) is 
identical to that of the GVB wave function, and 
the anti-bonding covalent spatial function |$^) 



from Equation (6a) possesses ungerade symme- 
try to couple to the m^ = triplet spinor. 
Thus, the variational parameter A, which be- 



comes non-zero for all bond lengths greater 
than the CF point, controls the spin symme- 
try breaking. In the limit where A ^ 1 these 
orbit als simply reduce to (pi and 0^ above; in the 
limit A ^ the wave function becomes identi- 
cal to RHF; and in the limit A ^ — 1 the wave 
function describes a spin-contaminated charge 
density wave. While the UHF wave function 
correctly dissociates the molecule using a single 
determinant, this comes at a rather high price: 
a discontinuity in the derivative of the energy 
with respect to the internuclear separation oc- 
curs at the CF point, and at large distances the 
wave function no longer possesses the correct 
parity or spin. 

It is, of course, possible to break further 
symmetries of a single-determinant wave func- 
tion. These have been explored in detail by 
Fukutome^^ as well as Stuber and Paldus.'^S 
In our opinion, this work has not received 
the attention that it deserves because break- 
ing additional symmetries sacrifices more quan- 
tum numbers that, in turn, become quite dif- 
ficult to recover. An approach of particular 
note, the generalized Hartree-Fock approxima- 
tion (GHF), allows each electron to occupy an 
orbital without definite a or /3 spin. Computa- 
tionally, one typically accomplishes this by al- 
lowing the MO coefficients to become complex, 
i.e. each real atomic orbital gets both a real 
and imaginary a and /3 coefficient. The GHF 
approach is fraught with other difficulties, as 
dissociation curves often contain energy cross- 
ings and other weird behavior. For these rea- 
sons, the community has largely neglected the 
GHF approach. However, preserving some sym- 
metries of the wave function in a single determi- 
nant approach (e.g. UHF) comes with an addi- 
tional price: closed shell molecules can't in gen- 
eral be dissociated into the correct open-shell 
fragments. However, using the GHF approach 
recovers size consistency in the resulting wave 
function.^ We direct the reader to the rather 
vast literature on the subject of geminal func- 
tions and symmetry for an in depth discussion 
(see, for example, ref. ( [20]) and ( |23j) and the 
references therein), and have simply presented 
a few conceptual ideas here which shed light on 
some of the issues that arise when we encounter 



spin projection below. 

2.3 Pairing in Broken Symmetry 
Wave Functions 

The UHF wave function possesses a particu- 
larly interesting property, pointed out in differ- 
ent forms by many authors. Lowdin proposed 
the "pairing" theorem ^^ which states that the 
orbitals of a wave function of the UHF type 
can always be transformed such that the alpha 
and beta orbitals are mutually orthogonal ex- 
cept for "corresponding pairs" of orbitals which 
have non-zero overlap. This follows from the 
corresponding orbital transformation of Amos 
and Hall^^ or, alternatively, from the left and 
right eigenvectors of a singular value decompo- 
sition. 

Consider a UHF wave function defined by a 
set of occupied alpha spin orbitals, ip^^ and a set 
of occupied beta spin orbitals, ip^ ^ orthonormal 
within themselves, but not mutually so. De- 
noting overlap between the two sets as S^f^ 
they can be brought to bi-orthogonal form by 
diagonalizing the matrices SS^ and S'^'S withing 
the subspace of the occupied a and /3 electrons, 
respectively, to generate two new sets of spin- 
orbitals, JJip"^ and V?/^^: 

SS^U = CU and S^SV = C'V. (18) 

The overlap matrix becomes diagonal in these 
subspaces, i.e.: 



UtSV = 5i,v^., 



(19) 



and the eigenvalues ^i = ^^ are bound between 
and 1. Thus, each a orbital pairs with a single 
/3 orbital, and the overlap for pair i is given by 

Although a somewhat straightforward result, 
it does not seem to have been emphasized that 
because the corresponding orbitals U and V 
are obtained from unitary transformations of 
the occupied a and /3 UHF orbitals, the UHF 
wave function may be written as an APSG: 
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ij = A[core]Y[uj,{l,2), 



(20) 



i=l 



where the geminals cl;^(1, 2) are given by: 



(21) 



and the spatial overlap of the non-orthogonal 
alpha and beta orbit als making up the pair is: 



bi 



{Ui\vj) = y/^i 



(22) 



This pair function extends the GVB pair idea, 
retaining the condition of strong orthogonality, 
but relaxing the perfect pairing constraint and 
allowing spin contamination from the triplet 
component of the state a/3 associated with each 
geminal. The non-orthogonal representation of 
the pair above can of course be transformed 
back into an orthogonal "natural orbital" rep- 
resentation, or expressed in the space of the 
symmetry restricted molecular orbitals as in 
Equation [8| Thus the UHF wave function also 
contains only intra-pair correlation, though re- 
laxing the perfect pairing constraint allows it 
the ability to capture more correlation than the 
GVB approach at the cost of the spin quan- 
tum numbers. Finally, we wish to emphasize 
that the GHF wave function is not an APG, as 
it also captures inter-pair correlation, retaining 
size consistency. 

2.4 The Broken Symmetry IPDM 

Lowdin pointed out the usefulness of the one- 
particle density matrix (IPDM) for defining the 
natural spin-orbitals of a system, and noted 
that the eigenvalues of this matrix were con- 
strained to be between zero and one.'^Sl Cole- 
man later argued that all two-particle interac- 
tions can be described using the one and two- 
particle density matrices. ^^ We therefore pause 
for a moment to remark on the structure of the 
IPDM of the UHF and GHF wave functions. 

One can separate the IPDM into spin blocks 
and write the IPDM of a UHF wave function 
in the following way: 
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(23) 
(24) 



The half sum and half difference of the non-zero 
blocks form the charge (P) and spin densities 
(M^), respectively: 



Paa + P(3^ 



and 



M, = 



P/3^ 



(25) 
(26) 



The magnetization describes the UHF spin con- 
tamination arising from the triplet contribution 
to the broken symmetry wave function. We call 
such a wave function collinear because Sz re- 
mains a good quantum number. 

The IPDM associated with the UHF wave 
function must be both Hermitian and idempo- 
tent since this wave function is a single Slater 
determinant. One can then trivially show the 
idempotency of each each spin block Pcrcr- Thus, 
the charge density is a linear combination of 
two idempotent matrices, and consequently has 
eigenvalues which are either 0, 1, 1/2, or come 
in "corresponding pairs" of n and 1 — n.^^^^The 
natural orbitals of the charge density coupled as 
"corresponding pairs" are simply the natural or- 
bital representation of the non-orthogonal UHF 
geminals discussed above. 

In the GHF formalism, even though addi- 
tional symmetries are broken, the wave func- 
tion is still a single Slater determinant though 
the IPDM takes the more general form: 



Paa Pa(3 
P/3a Pl3/3 



P + M^ M^ - iMy 



M, + M, 



M. 



(27) 
(28) 



where the x and y components of the magneti- 
zation are given by: 



M. = 



M„ 



Pa/3 + Pl3a 



Pal3 — Pl3o 



2i 



and (29) 
(30) 



For a GHF wave function, no orbital rotation 
can bring the IPDM to the form of Equation 



(23), and we call such a state non- collinear^ as 



the state no longer has a good quantum number 
associated with Sz- 
The most general GHF wave function can- 



not be written as the sum of two idempotent 
matrices, and loses the property of correspond- 
ing pairs. This results from the fact that the 
wave function cannot be written as an APG. 
Where the UHF formalism can only describe 
intra-pair correlation, the GHF formalism has 
additional flexibility to capture inter-pair cor- 
relation. Since the pairing in the IPDM is 
between alpha and beta spins, the UHF ap- 
proach only describes opposite-spin correlation, 
whereas through inter-pair coupling the GHF 
approach possesses the flexibility to describe 
same-spin correlation as well. 

2.5 Projected Hartree-Fock 

The primary drawback of broken symmetry 
wave functions lies in the difliculty of recap- 
turing quantum numbers associated with the 
broken symmetries. Many attempts have been 
made to address this problem ,ESESI originat- 
ing with the work of Lowdin.^ Recently a 
signiflcant breakthrough was made in deter- 
mining, self-consistently, symmetry adapted 
wave functions from a single broken-symmetry 
(deformed) determinant.'^ The resulting wave 
function possesses multi-reference character 
calculated with mean- field computational cost. 
This projected quasi-particle theory (PQT) 
relies upon the simple idea of deliberately 
breaking and then self-consistently restoring 
the symmetries of the Hamiltonian with the 
use of projection operators. A computationally 
eflBcient formulation of projected Hartree-Fock 
theory based upon PQT recently appeared in 
the literature. ^^ In that work, breaking and 
restoring discrete symmetries such as point 
group and complex-conjugation symmetries 
were also addressed. In this work, we focus on 
spin symmetries as they are most familiar to 
quantum chemists, though the method employs 
an approach analogous to that of restoration of 
angular momentum in the context of nuclear 
physics. ^^ 

Let us begin by taking a moment to review 
the two types of wave functions we will project. 
The collinear UHF wave function \UHF) min- 



imizes the energy functional: 



E 



(*l^l*) 
(^1^) ^ 



(31) 



sacrifices spatial and spin symmetry, and re- 
tains Sz as a good quantum number. The 
non-collinear GHF wave function minimizes the 
same functional with the additional sacrifice of 
the Sz quantum number. We will denote pro- 
jections of these two wave functions as \SUHF) 
and \SGHF)^ respectively. We cannot empha- 
size enough that these projected wave func- 
tions are not single Slater determinants. Con- 
sequently, the density matrices associated with 
these multi-reference wave functions are not 
idempotent, though they do preserve the sym- 
metries of the Hamiltonian. 

To determine our spin-projected wave func- 
tions, we write the energy of a projected, de- 
formed determinant 1$) as: 



E 



($|P|$) 



(32) 



where H is the Hamiltonian, and P is the 
(Hermitian and idempotent) projection opera- 
tor that restores the symmetries broken in the 
deformed determinant |$). Formally, for con- 
tinuous symmetries such as S'^ and Sz^ the pro- 
jection operator takes the form of an integral. 
In practice, we discretize this integral over a 
grid of modest size. 
It can then be shown that the energy expres- 



sion (32) can be written as a functional of the 



one-particle density matrix, p<i>, associated with 
the underlying deformed determinant |$). One 
then minimizes the resulting functional with 
respect to the deformed orbitals. The result- 
ing PHF equations are qualitatively similar to 
that of Hartree-Fock and the method retains 
mean- field computation cost. Further, as with 
Hartree-Fock, any symmetry present in the ini- 
tial guess for the deformed determinant will be 
preserved throughout the optimization proce- 
dure .l^^l 

We pause for a moment to remark that the 
implementation described in ref. ( 24) does not 
strictly restore aS^, as this would require a two- 



body operator. Instead, it restores rotational 
invariance of the deformed determinant in spin- 
space. For coUinear determinants, we rotate the 
determinant so that it hes along the z-direction, 
and then make it rotationally invariant using a 
projection operator of the form: 
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(33) 



where (i^^(/3) = {s]m\e'^^^y\s]m). This leaves a 
particular direction as "special" in that the pro- 
jected wave function (and its energy) depends 
on the choice of the quantum number m = {Sz) 
of the coUinear determinant. Projecting a non- 
collinear determinant requires integration over 
three Euler angles and eliminates the depen- 
dence on the quantum numbers of the deformed 
state. 

Interestingly, the projected wave function, 
though multi-reference in nature, may be char- 
acterized by a single, deformed determinant we 
denote as |$) or, equivalently, the associated 
one-particle density matrix p^. We want to 
stress that |$) ^ \UHF) or \GHF). The de- 
terminants \UHF) or \GHF) optimize the en- 
ergy functional of Equation ( [31] ) (breaking ei- 
ther S'^, or S'^ and S^^, respectively), whereas 
the determinant |$) optimizes the energy func- 
tional of Equation (32). The determinant |$) 



can be of either the UHF or GHF type, depend- 
ing on what symmetries are broken. 

Finally, we remark that though the SUHF 
wave function and the spin projected extended 
Hartree-Fock (EHF) method are in fact equiv- 
alent, we obtain our results without working 
in the corresponding orbital basis (as suggested 
by Mayer and Lowdin) resulting in significantly 
less computational overhead.!^ In addition, pre- 
vious examples of PHF in the literature have 
all been coUinear, i.e. restoring only S^ ^ rather 
than both S^ and Sz^ i.e. previous work did 
not consider SGHF wave functions. As Pu- 
lay pointed out,!^ when several orbitals are 
strongly correlated, the EHF (equivalently the 
SUHF) wave function doesn't have the fiexibil- 
ity to describe all of the static correlation in 
the system. The SGHF wave function provides 
a way around this: by breaking and restoring Sz 



as well as aS^, it captures inter-pair correlation 
which the SUHF wave function cannot. 



3 Hydrogen Networks 

3.1 Computational Details 

We utilize an implementation of the PHF equa- 
tions in the development version of the GAUS- 
SIAN suite of programs. ^^ A minimal (ST0-3G) 
basis set was deemed sufficient for the purpose 
of describing the qualitative features of strong 
correlation associated with near degeneracies in 
the orbital spectrum. The initial guess for de- 
formed SUHF determinants are broken symme- 
try Hartree-Fock states. The initial guess for 
deformed SGHF calculations were constructed 
using an SUHF determinant and mixing the al- 
pha and beta orbitals nearest the Fermi level 
with a small mixing angle. 



3.2 H4 and H^ Rings 

We illustrate the evolution of pair wave func- 
tions as a function of nearest-neighbor distance 
{Rn-n)-> as well as investigate the errors asso- 
ciated with different types of spin projection 
(SUHF and SGHF), by considering the sym- 
metric stretching of rings of four and six hy- 
drogen atoms. There is a topological difference 
between the orbital spectra of the four and six 
membered rings, and in general between 4N and 
4N+2 membered rings, indicated in Figure [2] 
In general, 4N membered (anti-resonant) rings 
have an effective Fermi level which lies within 
two half-filled degenerate orbitals, whereas the 
effective Fermi level for the 4N+2 membered 
(resonant) rings lies between a pair of fully oc- 
cupied degenerate orbitals and a pair of unoccu- 
pied degenerate orbitals. One expects that the 
ground state for both systems should be a sin- 
glet, though in the case of anti-resonant rings, 
an open-shell singlet. 

The ground state of H4 adopts a la^(ea;,e^)^ 
configuration, with x and y defined in the in- 



set of Figure 3a, The open-shell configuration 
gives rise to three singlet states: ^Ag{x^ + 2/^), 



ngy^ —y^)i and a triplet state: 
^B2g{xy). Placing the electrons in ^Big{x'^ ^'^^ 



^B2g{xy), and ^B^n[x^ 



-y') 



--Ef- 



H4 



He 



Figure 2: A schematic depiction of the one- 
electron molecular orbital energy spectrum for 
the H4 (left) and He (right) rings, with Lowdin's 
pairing scheme indicated by arrows. The effec- 
tive Fermi energy (E^?) lies in the degenerate 
orbitals for H4 and in the gap in Hq. 

yields the lowest energy. The reason can be 
understood in the following way: the ^Big 
state takes the form of an open-shell singlet 
^l+y^ e^_^, where one electron resides on each of 
the atoms lying on one diagonal of the square 
(x+y), and the other along its complement (x- 
y). This segregates the alpha and beta spins 
onto separate sublattices allowing them to avoid 
double occupancy and the consequent on-site 
repulsion. The electrons in the lowest (1%) 
orbital then correlate by mixing with its anti- 
bonding counterpart {Ibg) as shown in Figure 

El 

We present the ^Big dissociation curves 
for several approximations in Figure [3a| as a 
function of nearest-neighbor distance (Rn-n)- 
While the RHF approximation dissociates im- 
properly, both the UHF and GVB approxima- 
tions succeed. Figure [3b] shows the energy dif- 
ference from the exact (FCI) result and one can 
see that, in the intermediate coupling regime 
{Rn-n < 2.0A), UHF has a significantly smaller 
error than GVB. This refiects the advantage of 
allowing orbital spin contamination in the UHF 
approach. Neither of these wave functions are 
equal to the FCI for this system, and the error 
in the absolute energy is on the order of 1-2 eV 
near equilibrium. 

In the intermediate recoupling regime where 
both UHF and GVB still have significant er- 
ror, the projected wave functions offer signifi- 
cant improvement in the energy, in addition to 
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Figure 3: (a) Dissociation curves for different 
singlet wave functions for the ground state of 
the H4 square in a minimal basis. The zero 
of energy has been set at the energy of four 
Hartree-Fock hydrogen atoms, (inset) Arrange- 
ment of the atoms with respect to x and y sym- 
metry axes. |(b)| Energy difference from FCI 
of several singlet wave functions for H4. The 
SUHF energy remains within 3 meV of the ex- 
act solution at all distances. 



restoring the correct quantum numbers. This 
improvement results from the multi-reference 
nature of the projected wave function. The 
error in the SUHF wave function reduces to 
the order of meV, and the SGHF wave func- 
tion is identical to the FCI for the minimal ba- 
sis. Recall that for H2 both the SUHF and 
the GVB-PP approximations are identical to 
the exact wave function for the minimal basis 
set. While they continue to refiect the appro- 
priate intra-pair correlation in the four electron 
system, they do not treat the inter-pair corre- 
lation completely, a problem ultimately associ- 
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ated with the lack of size-consistency in these 
methods. 

To understand these results, recall that there 
exist two distinct singlet spin states for any 
four-particle fermionic systemiEU 



|15) = |0,0)|0,0) 

1 
1 



\2S) 



1,+>|1, 



+ l^)|l,-)|l,+) 



(34a) 



(34b) 



^)|1,0)|1,0). 



Here, \1S) describes the product of two singlet 
coupled pairs as in the GVB-PP wave func- 
tion, while \2S) describes a linear combination 
of products of triplet pairs coupling to form 
an overall singlet state. The dominant contri- 
bution to the ground state comes from \IS)^ 
though the exact wave function does contain a 
small contribution from |25'), particularly near 
the equilibrium bond length. The GVB-PP 
wave function successfully describes all of the 
contribution to the total energy associated with 
the \1S) spin state. The UHF singlet obtains a 
better total energy by breaking spatial symme- 
try and introducting spin contamination. 

For the spin projected states, the SUHF 
wave function captures correlation from deter- 
minants that couple to the \1S) spinor. Though 
both GVB-PP and SUHF wave functions are 
multi-reference and both couple to |1aS), SUHF 
captures significantly more correlation because 
it contains more determinants than the GVB- 
PP. The first two terms of \2S) describe same- 
spin correlation, i.e. contributions from the 
product of m^ = ±1 two-electron triplets. 
There are no same-spin correlations present in 
a deformed rris = singlet determinant; con- 
sequently, the SUHF state cannot describe cor- 
relation from determinants that couple to |2S'), 
resulting in the energy difference from FCI seen 
in Figure |3b| Breaking Sz as well as S'^^ allows 
the SGHF wave function to further correlate 
the two pairs by capturing this same-spin cor- 
relation. By deliberately breaking and restoring 
this symmetry, the SGHF wave function recov- 



ers the full inter-pair correlation energy, i.e. the 
FCI result. 



0.3 



0.2 



0.1 



0.0 



-0.1 



UHF 

GVB 3g 

SUHF m, = 1 

SUHF m = 



2g 






1.0 



1.5 



2.0 

Rn-n(A) 



2.5 



3.0 



(a) 



16.0 
14.0 
12.0 - 
10.0 - 

8.0 

6.0 

4.0 - 

2.0 - 

0.0 



"Mo" 

60.0 

30.0 

0.0 



----'* 



1.0 1.5 2.0 2.5 3.0 

.....•••■■'■' UHF 

..••••■'■" GVB 3p. 

SUHFm,= l --■ 
::;: SUHF m = 



1.0 



1.5 



2.0 

Rn-n(A) 



2.5 



3.0 



(b) 



Figure 4: (a) Dissociation curves for an H4 



square for different triplet wave functions in a 
minimal basis. The zero of energy has been 
set at the energy of four Hartree-Fock hydrogen 
atoms. |(b)| Energy difference from FCI for the 
same wave functions. 

There is another point we wish to make with 
this simple model exemplified in the triplet 
potential curves depicted in Figure |4} The 
GVB-PP ^B2g{xy) triplet state arising from the 
ground state configuration has the wrong sym- 
metry to dissociate properly as evident in Fig- 
ure [4al Another triplet state, based on the ex- 
cited electronic configuration la^(ea;,e^)^ (not 
shown) does, but the dissociation curve for this 
state is purely repulsive and not appropriate at 
equilibrium. The UHF solution, by breaking 
spatial symmetry, yields a qualitatively correct 
curve at all distances, though it predicts a bar- 
rier to formation of the bound state. 

Turning now to the projected wave functions. 
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we see that, as in the singlet case, the SGHF 
triplet wave function is equal to FCI for this 
minimal basis. However, we find two distinct 
SUHF triplet wave functions with very differ- 
ent energies; one arises from projecting a de- 
formed determinant with m^ = and a second 
arises from projecting a deformed determinant 
with rris = ±1. As we briefiy mentioned above, 
these two states differdue to the fact that we 
simply assure rotational invariance of the pro- 
jected state, without respect for a particular 
axis in spin space. When we restore both ro- 
tational invariance as well as Sz (i.e. SGHF) 
we don't have this problem. 

Using this simple four-electron system, we 
can understand why this occurs. The PHF 
equations preserve the symmetries of the ini- 
tial guess throughout the SCF procedure. The 
symmetries of the two deformed determinants 
are diflFerent and thus capture different types of 
correlation. Recall that there exist three dis- 
tinct four-particle triplet spinors:'^ 



I IT) 
\2T) 

|3T) 



|0,0)|1,0) 
|1,0)|0,0) 



(35a) 
(35b) 

|1,-)|1,+)) (35c) 



We obtain the SUHF m^ = 1 wave function 
by taking the UHF open-shell triplet state as 
an initial guess. For distances smaller than 
the triplet CF point, the deformed determi- 
nant has a doubly occupied lag orbital and the 
open-shell electrons occupy the Cx and Cy or- 
bit als, delocalizing the un-paired spins across 
all the atoms. Beyond the triplet CF point, the 
deformed determinant breaks symmetry, mix- 
ing the Cx and Cy orbitals as described above 
to localize the open-shell spins on different 
atoms. We therefore see that beyond the CF 
point, the two pairs must form the product 
spinor |3T), as the open-shell geminal's spinor 
takes the form of either |1,+) or |1,— ), and 
the remaining geminal couples to the spinor: 
4= (|1, +) — |1, — )). This provides only a mod- 
est improvement in the energy compared to 
UHF as symmetry breaking captures most of 
the correlation. 



function state by taking the UHF singlet state 
as an initial guess. This deformed determi- 
nant captures contributions from |1T) and |2T), 
and the resulting SUHF state captures determi- 
nants which couple to product geminals from 
the open-shell singlet and the paired triplet 
(and vice versa). As there are significantly more 
determinants in the FCI space which couple to 
these four-particle spinors, the maximum error 
drops from almost 2 eV to slightly less than 90 
meV (see Figure 4b, inset). 
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In contrast, we obtain the SUHF m^ = wave 



Figure 5: The singlet-triplet energy difference 

for H4. Sharp cusps in the UHF and SUHF 

rris = I states are the result of spontaneous 

symmetry breaking at a CF point. 

This has consequences for calculating mag- 
netic properties of materials. We plot the 
singlet-triplet splitting for the H4 system in Fig- 
ure [5} At large Rn-n this splitting relates to the 
magnetic coupling constants of the Heisenberg 
model Hamiltonian. The dissociation curves of 
both the UHF and the SUHF rris = I state have 
a CF point, and consequently the singlet-triplet 
splitting calculated from these states possesses 
a discontinuity in the derivative at that point. 
Combined with the significant underestimation 
of the triplet energy, the splitting calculated us- 
ing the SUHF rris = I wave function is actually 
worse than using the UHF single determinant. 
In contrast, the rris = SUHF state gives both 
a quantitatively as well as qualitatively correct 
curve. 

Turning now to Hq^ two additional electrons 
fill the degenerate e^ orbitals, leading to a 
closed-shell singlet ground state. In Figure [6| 
we plot the UHF pair functions for the three 
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Figure 6: UHF corresponding (pair) orbit als for the He ring at short, intermediate and large 
separation. ^^ We present the pairs with the a orbital on the left and the corresponding /3 orbital 

on the right; ^/^ is the orbital overlap (see text). 



pairs of Hq ring at three distances: 1.2 A, 1.8 
A, and 5.00 A, corresponding to just past the 
CF point, moderately stretched, and dissocia- 
tion. The pairs evolve smoothly as the a and 
P electrons segregate onto two distinct sublat- 
tices. In order to understand the manner in 
which the pair functions evolve, we project the 
pair functions back onto the symmetry adapted 
MOs of the system. Orbital coefficients for the 
alpha and beta pairs for this system^^l in the 
symmetry adapted MO basis at a distance of 
1.80 A appear in Table [l} showing how the 
pairs form by mixing a bonding orbital with 
its anti-bonding counterpart. In the dissocia- 
tion limit, the orbital coefficients approach the 
limiting value of 1/a/2- The pairing ansatz ev- 
ident in the UHF wavefunction, where each or- 
bital mixes with its corresponding anti-bonding 
orbital, was suggested as a way to correlate elec- 
trons subject to large Coulomb interactions by 
Lowdin in 1955 and goes by the name of alter- 
nant molecular orbital (AMO) theory."^^^ For 
the projected wave functions, the underlying 
deformed determinant of the SUHF state pos- 
sesses pairing scheme. 

To summarize, we believe these approaches 
capture the most important static correlation 
effects. With the SUHF wave function, we cap- 
ture all the intra-pair correlation effects, and 



the SGHF approach captures additional pair- 
pair correlation. Even so, higher order many- 
body effects are not described properly by these 
methods, and those terms grow in importance 
as the system size increases, leading to a lack 
of size extensivity.^ 

3.3 Larger Rings 

We now consider larger ID rings of hydrogen as 
a function of nearest-neighbor bond length. We 
present results for both 16 and 18 atom rings; 
chosen because they are complicated enough 
that FCI calculations become quite cumber- 
some, yet the system remains small enough that 
the PHF approach can provide a significant im- 
provement over UHF, i.e. we do not approach 
the thermodynamic (A^e -^ co) limit. Our re- 
sults simply extend those already seen in the 
four and six atom rings, and in what follows we 
focus on the resonant system: H^^g. 



Figure 7a shows the UHF one-electron molec- 
ular orbital energy spectrum for the H^i^g ring. 
In the small Rn-n regime, all the orbit als except 
for the lowest occupied and highest virtual or- 
bitals remain doubly degenerate. As the bonds 
stretch, the orbital spectrum compresses as the 
occupied and virutal MOs become energetically 
degenerate. At dissociation, the orbital eigen- 
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Table 1: Corresponding pair coefficients in the symmetry adapted basis'^ 
(1.80 A) bond length. Orbitals pair as shown in Figure 
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values collapse into two sets of degenerate lev- 
els, an occupied "lower band" and the unoccu- 
pied "upper band", each containing nine spatial 
orbitals. The effective Fermi energy sits in the 
gap between these two "bands." 

Although our calculation uses the exact 
Hamiltonian, we map our results as a function 
of distance to a Hubbard model in an attempt 
to make contact with the substantial body of 
work using that approximation. To this end, we 
approximate the onsite interaction parameter U 
as the UHF "band gap" in the limit Rn-n -^ co. 
In the case of half-filling, this energy difference 
corresponds to the energy needed to move an 
electron from one site to another. This "gap" 
has a value of ^ 21 eV for both the H^^g and 
H^g rings, not an unreasonably large value for 
U in this minimal basis set. 

We next extract the effective near-neighbor 
hopping parameter t from the band width of 
the one-electron spectrum. The effective value 
we seek should be renormalized by all the terms 
in the exact Hamiltonian except for the on-site 
U term. To this end, at each bond length we 
compute the UHF eigenvalue spectrum for the 
high-spin, ferromagnetic state. This singly oc- 
cupies each orbital, thereby eliminating the on- 
site Hubbard U interaction, and the effective 
hopping parameter t is then given by: 



* ^ ^ (^A^ - ^o) 



(36) 



where e^ - cq is the splitting between the high- 
est and lowest occupied MOS of the ferromag- 
netic state, and z is the number of nearest 
neighbors. 



The onset of strong correlations corresponds 
to the regime where the bandwidth and the 
on-site repulsion become comparable. For the 
ID ring, this occurs at a ratio of U/t ^ 4, 
or Rn-n ^ 1-25 A, shortly past the CF point. 
As Rn-n -^ CO the effective hopping parameter 
tends toward zero, refiecting an effectively infi- 
nite on-site interaction. This transition regime 
corresponds to distances where the SUHF en- 
ergy differs most greatly from the exact result 



for H4 (the inset of Figure 3b) suggesting the 
importance of inter-pair correlation. 
In the dissociation curves for the H^g ring 



shown in Figure [7b| the UHF wave function 
goes to the correct limit by allowing the al- 
pha and beta spins to localize on different sub- 
lattices (See Fig 7b, inset). Notice that the 
CF point is quite near equilibrium, i.e. the 
UHF wave function starts to correlate through 
localization in this region. One might there- 
fore suspect that strong correlations contribute 
significantly at these distances. In fact, the 
SUHF projected wave function captures ^3.5 
eV of energy from intra-pair correlation, and the 
SGHF wave function captures an additional ^ 
1.4 eV from inter-pair correlation at the equi- 
librium near-neighbor distance. 

The corresponding pairs in H^g and H^g are 
similar to those discussed for H4 and He, respec- 
tively. At small distances, the orbital overlap 
of each pair remains identically one (see Fig- 
ure [s]) . This follows from the fact that inside 
the CF point, the UHF wavefunction has not 
yet broken symmetry and each spin-orbital {a 
and /3) remain spatially identical. The electrons 
nearest the Fermi level localize first as the ring 
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Figure 7: 



a) The UHF orbital energy spectrum 

and U/t. 



for the H^i^g ring as a function of R^ 
As the bonds are stretched, the occupied and 
virtual MOs each become energetically degen- 
erate, forming upper (unoccupied) and lower 
(occupied) Hubbard "bands." |(b)| Dissociation 
curves and (inset) UHF spin density for the H^g 
ring. The UHF wave function localizes the elec- 
trons anti-ferromagnetically at large Rn-n- 

stretches, and at dissociation the overlap of all 
the orbitals goes to zero, reflecting compete lo- 
calization of the spins onto separate sublattices. 
Again, we pause to stress that the multi- 
reference SUHF wave function can be charac- 
terized by a single, deformed determinant, dis- 
tinct from the self-consistent UHF wave func- 
tion. To illustrate this point, we show the 
value of the overlap integral {^/i) for core and 
degenerate HOMO orbitals for both the UHF 
and broken-symmetry determinant character- 
izing the SUHF wave function in the inset of 
Figure [8| By deliberately breaking the symme- 
try of the underlying determinant in the SUHF 
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Figure 8: Overlap of the singlet UHF corre- 
sponding orbitals for the H^^s ™S (see text). All 
orbitals except for the core are doubly degen- 
erate. The UHF overlap becomes identically 1 
for all orbitals at the CF point, (inset) Overlap 
of the a and /3 HOMO and core corresponding 
orbitals of the UHF and deformed SUHF deter- 
minants of the same system. 

approach, artiflcially localizing the electrons at 
small distances, and projecting the deformed 
determinant, the SUHF approach captures a 
signiflcant amount of correlation energy near 
equilibrium. 

Finally returning to the anti-resonant H^g 
ring, we plot the core and HOMO pair orbitals 
in Figure JQ) These are analogous to the case of 
H4, and the HOMO orbital localizes one of each 
spin on separate sublattices at all distances. 
This creates an open-shell singlet, while the rest 
of the orbitals retain the freedom to localize or 
not as correlation effects demand. Near the CF 
point the core orbital remains delocalized, and 
by 3 A it localizes each spin onto different sub- 
lattices. The remaining orbitals behave in a 
similar manner. 

3.4 2D Hydrogen Networks 

We conclude with a brief discussion of a 2D sys- 
tem: symmetric dissociation of H^g in a square 
geometry. The square lattice has been previ- 
ously studied at half-fllling by Larson and Thor- 
son, and Calais and coworkers. ^^^^ Though 
the system studied here doesn't have periodic 
boundary conditions, it still captures the es- 
sential Hubbard physics. We point out that 
the UHF single particle energy spectrum differs 
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(a) core pair at 1 A 



(b) core pair at 3 A 
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(c) HOMO pair at 1 A 



(d) HOMO pair at 3 A 



Figure 9: Core and HOMO UHF corresponding (pair) orbit als for the YL^q ring. At large Rn-n the 
broken symmetry wave function locahzes the alpha and beta spins on different sublattices. 



from that of the H^i^g ring in that the degeneracy 
of the HOMO has lifted and the Fermi energy 
now lies within the gap. This results in a close- 
shell singlet ground state, which becomes the 
focus of our discussion. 
Interestingly, we find two distinct stable UHF 



states, shown in Figure 10a The difference be- 
tween these two states is subtle, amounting to 
a ninety degree rotation of the localized spins 



in the center of the plaquette (see Figures 10b 



and 10c). We also plot the pair orbitals for the 
AFM state near dissociation,'^ showing the seg- 
regation of respective spins onto separate sub- 
lattices. The UHF global minimum energy oc- 
curs for the AFM state, however, the orbital 
Hessian indicates that the SF state exists as a 
true local minimum, and not a transition state. 
One might naively hope that projecting these 
two different UHF singlet states would yield the 
same multi-reference wave function for the sin- 
glet state. Unfortunately, this intuition fails for 
the SUHF approach. 

As mentioned previously, the SCF procedure 
preserves the symmetries of the initial guess 
for the broken symmetry determinant; conse- 
quently, we find two distinct SUHF singlets. 
Using the SGHF approach overcomes this prob- 
lem: breaking Sz allows the dissociation curve 
to smoothly interpolate between these two dif- 
ferent spin densities. The SGHF wave func- 
tion captures the most correlation in the region 



where these two states become energetically de- 
generate, as one would expect. Note that we 
have not performed a projection to restore spa- 
tial symmetry; the SGHF wave function only 
restores S'^ and S^- 

We also find multiple UHF and SUHF so- 
lutions arising from broken spatial symmetries 
in larger models such as 8x8 plaquettes. One 
can in some sense think of these different states 
as incorporating different types of spin fiuctu- 
ations. By interpolating between these spin 
states, the SGHF approach captures additional 
correlation in the region where states of differ- 
ent spatial symmetry become energetically de- 
generate. Ongoing work indicates that break- 
ing and restoring point group symmetry often 
becomes very important in cases such as these 
where the SGHF wave function does not have 
enough fiexibility to describe all of the static 
correlation in the system. 

We conclude this section by pointing out that 
we may apply the pairing ansatz which seg- 
regates spins onto separate sublattices to any 
alternant lattice: the one-dimensional, square, 
simple cubic and bcc lattices. For example: 
in terms of the Bloch functions of the orig- 
inal 2D lattice, the pairs mix the configura- 
tions iplx,ky with their anti-bonding partner 
^kx+7Tky+7T' The wave function may be writ- 
ten as a spin-contaminated APSG if the corre- 
sponding pairs at each k-value come from the 
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Figure 10: The H^^e plaquette. (a) Dissociation curves for singlet wave functio ns (b) AFM state 
UHF spin density, and (c) SF state UHF spin density at Rn-n = 3.0 A. (d) (k) UHF pair orbitals 
and overlaps for the singlet AFM state at Rn-n = 3.0 A.S21 At this distance, each pair localizes one 

alpha and one beta spin on different sublattices. 



corresponding orbitals of the spin-polarized so- 
lution. The energy of this wave function is iden- 
tical with that of the Neel state. In fact, the 
energies of all the pure spin states are degener- 
ate, and the Neel state a linear combination of 
them. 

4 Conclusion 

We have emphasized that the UHF wave func- 
tion can be written as a spin-contaminated pair 
wave function of the APSG form. The overlap 
of the alpha and beta corresponding orbitals 
of the UHF solution can be taken as a proxy 
for the strength of the correlation captured by 
breaking symmetry. As a function of distance, 
or the ratio U/t, the UHF corresponding or- 
bitals pair in a manner allowing a smooth evo- 
lution from a system with doubly occupied or- 
bitals into one in which the a and /3 electrons 
segregate onto distinct sublattices. In this way 
the UHF wave function evolves from a regime 



in which the material behaves as a metal with 
delocalized spins, to one in which all the spins 
are localized antiferromagnetically. 

Projecting the spin-contaminated UHF pair 
wave function recovers additional correlation 
energy in finite systems. The resulting multi- 
reference wave function, characterized by a sin- 
gle, deformed determinant, describes a pure 
spin state. The SUHF wave function (project- 
ing aS^ by assuring rotational invariance in spin- 
space) captures additional intra-pair correlation 
beyond UHF. The SGHF wave function (pro- 
jecting S'^ and Sz) captures additional inter- 
pair correlation. By deliberately breaking and 
then restoring these symmetries, we describe 
strong correlation even at geometries where the 
UHF solution reduces to RHF or the GHF so- 
lution reduces to UHF. 

From simple calculations of one and two- 
dimensional finite lattices of hydrogen atoms, 
we have shown that the broken symmetry de- 
terminant characterizing our SUHF solutions 
mirror the pairing schemes utilized extensively 
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some years ago in AMO theory to describe cor- 
relation in alternant pi-bonded networks. This 
early work did not address the questions of 
what happens as you dope the system away 
from half-filling, or in what manner the pairs 
respond. The pairs described here accounting 
for magnetic correlations are strongly orthogo- 
nal to one another, whereas the pairs of BCS 
theory overlap; nevertheless, the proximity of 
superconducting behavior to antiferromagnetic 
phases suggest this may be an interesting av- 
enue to explore. 
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